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ABSTRACT 

We reconsider the problem of gravitational structure formation inside and outside General Relativity 
(GR), both in the weakly and strongly non-linear regime. We show how these regimes can be explored 
observationally through clustering of high order cumulants and through the epoch of formation, abun- 
dance and clustering of collapse structures, using Press-Schechter formalism and its extensions. We 
address the question of how different are these predictions when using a non-standard theory of Gravity. 
We study examples of cosmologies that do not necessarily obey Einstein's field equations: scalar-tensor 
theories (STT), such as Brans-Dicke (BD), parametrized with to, a non-standard parametrisation of the 
Hubble law, H 2 = a~ 3 ( 1+e \ or a non-standard cosmic equation of state p — 7/5, where 7 can be chosen 
irrespective of the cosmological parameters (Q,m and ^2a)- We present some preliminary bounds on 7, lu 
and e from observations of the skewness and kurtosis in the APM Galaxy Survey. This test is indepen- 
dent of the overall normalization of rms fluctuations. We also show how abundances and formation times 
change under these assumptions. Upcoming data on non-linear growth will place strong constraints on 
such variations from the standard paradigm. 

Subject headings: Large-scale structure of Universe - galaxies: formation - Gravitation - instabilities 



1. INTRODUCTION 

In Cosmology the standard picture of gravitational 
growth, and also many aspects of fundamental physics, 
are extrapolated many orders of magnitude, from the 
scales and times where our current theory of gravity (Gen- 
eral Relativity, GR) has been experimentally tested, into 
the distant universe. In particular, current limits on the 
(parametrized) Post Newtonian formalism mostly restrict 
to our very local Universe (see Will 1993). It is important 
to evaluate how much our predictions and cosmological 
picture depend on the underlying hypothesis (see Peebles 
1999 for insightful comments on the state of this subject). 
The other side of this argument is that cosmology can be 
used to test fundamental physics, such as our theory of 
gravity. 

One aspect of GR that could be questioned or tested 
without modifying the basic structure or symmetry of the 
theory are Einstein's field equations, relating the energy 
content (T M j,) to the curvature (R^). One such modifi- 
cation, which will be considered here, is scalar-tensor the- 
ories (STT), such as Brans-Dicke (BD) theory. A more 
generic, but also more vague, way of testing the impor- 
tance of Einstein's field equations is to model indepen- 
dently the geometry and the matter content, thus allow- 
ing for the possibility of other relations between them. 
Some simple aspects of this idea will be illustrated here by 
studying structure formation in a flat, matter dominated 
universe but with a m ore g eneral growth law for the Hub- 
ble rate — see section [3.2| below. Similarly, we will also 
consider results for a generic equation of state: p — -fp, 
where 7 can be chosen independently of the cosmological 
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parameters (f^M, ^fc and 

Our aim in this paper is to explore certain variations of 
the standard model to see how they affect structure forma- 
tion. The idea is to find a way to parameterize variations 
from GR that might produce differences large enough to 
be observable. The variations considered could have other 
observable consequences (eg in the local universe or in the 
radiation dominated regime) which might rule them out 
as a viable new theory. But even if this were the case, we 
still would have learn something about how structure for- 
mation depends on the underlying theory of Gravity or the 
assumptions about the equation of state. This aspect of 
the theory has hardly been explored and it therefore repre- 
sents an important step forward in analyzing alternatives 
to the current paradigm, eg non-baryonic matter (see Pee- 
bles 1999), and could also help to set limits on variations 
of GR or the equation of state at high red-shifts. 

Here we consider two main regimes for structure for- 
mation in non-standard gravity/cosmology: weakly non- 
linear and strongly non-linear large scale clustering. We 
study the shear-free or spherical collapse (SC) model, 
which corresponds to the spherically symmetric (or local) 
dynamics (see below). This approximation works very well 
at least in two different contexts, that will be explored 
here. 

The first one is the growth of the smoothed 1-point 
cumulants of the probability distribution for large scale 
density fluctuations: the SC model turns out to repro- 
duce exactly the leading order perturbation theory pre- 
dictions (Bernardeau 1992), and turns out to be an excel- 
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lent approximation for the exact dynamics as compared 
to N-body simulations both with Gaussian (Fosalba & 
Gaztahaga 1998a, 1998b) and non-Gaussian initial con- 
ditions (Gaztanaga & Fosalba 1998). The measured 1- 
point cumulants in galaxy catalogues have been compared 
with these predictions (eg Bouchet et al. 1993, Gaztanaga 
1992,1994, Gaztahaga & Frieman 1994, Baugh, Gaztahaga 
& Efstathiou 1995, Gaztahaga 1995, Baugh & Gaztahaga 
1996, Colombi etal 1997, Hui & Gaztahaga 1999). 

The second one is the study of the epoch of formation 
and abundance of structures (such as galaxies and clus- 
ters), using the Press & Schechter (1974) formalism and its 
extensions (eg Bond et al. 1991, Lacey & Cole 1993, Shcth 
& Lemson 1999, Scoccimarro et al. 2000). Given some 
Gaussian initial conditions, this formalism can predict the 
number of structures (halos) of a given mass that will form 
at each stage of the evolution. One can use the SC model 
to predict the value of the critical linear over-density, <5 C , 
that will collapse into virialized halos. It turns out that 
the analytical predictions for the halo mass function and 
formation rates are remarkably accurate as compared to 
N-body simulations (Lacey & Cole 1994). One can also use 
this type of modeling to predict clustering properties of ha- 
los (eg Mo & White 1996, Mo, Jing & White 1997), cluster 
abundances (White, Efstathiou & Frenk 1993, Bahcall & 
Fan 1998) or weak lensing through mass functions (Jain 
& Van Waerbeke 2000). The observed cluster abundances 
have been used as a strong discriminant for cosmologi- 
cal models and also as a way to measure the amplitude 
of mass fluctuations, a% (see White, Efstathiou & Frenk 
1993, Bahcall & Fan 1998). 

In summary, we propose to address a very specific ques- 
tion here: how different are the above non-linear pre- 
dictions when using a non-standard cosmology and non- 
standard theory of Gravity? To answer this question we 
will consider two non-standard variations: scalar-tensor 
models and some examples of a cosmology that do not 
obey Einstein's field equations. The paper is organized as 
follows: In §2 we give a summary on how non-linear struc- 
ture formation relates to the underlying theory of Gravity 
(see Weinberg 1972, Peebles 1993, Ellis 1999 and refer- 
ences therein, for a review on the relation between grav- 
itational theory and cosmology). This section covers old 
ground with some detail as an introduction to later sec- 
tions and for the reader that is not familiar with this sub- 
ject or notation. We also present the more general case 
of an ideal (relativistic) fluid. As far as we know, some 
of the non-linear results presented here are new. In §3 
we show how these predictions change in the two exam- 
ples of non-standard gravity. Observational consequences 
are explored in §4. In §5 we present a discussion and the 
conclusions. 

2. GRAVITATIONAL GROWTH INSIDE GR 

The self-gravity of an over-dense region work against the 
expansion of the universe so that this region will expand at 
a slower rate that the background. This increases the den- 
sity contrast so that eventually the region collapses. The 
details of this collapse depends on the initial density pro- 
file. Here we will focus in the spherically symmetric case. 
We will revise non-linear structure growth in the context 
of the fluid limit and the shear-free approximation. These 
turns out to be very good approximation for the appli- 



cations that will be considered later (leading order and 
strongly non-linear statistics). 

We start with the Raychaudhuri's equation, which is 
valid for an arbitrary Ricci tensor R^ v . We use Einstein's 
field equations and the continuity equation to turn Ray- 
chaudhuri's equation into a second order differential equa- 
tion for the density contrast. We first present the matter 
dominated (non-relativistic) case, with sol utio ns for the 
linear and non-linear regimes. Later, in § |2.5| , we assess 
the more generic case of an ideal (relativistic) fluid and its 
corresponding solution. 

2.1. Einstein's and Raychaudhuri's Equations 

We start recalling that the metric tensor <? M „ defines the 
line element of space-time: 

ds 2 = dx» dx u (1) 

which in the homogeneous and isotropic model of the cos- 
mological principle can be written as (see eg Weinberg 
1972): 



ds 2 = dt 2 -a 2 {t) 



dr 2 



1 + kr 2 
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sin 



(2) 



As usual we will work in comoving coordinates x re- 
lated to physical coordinates by r p = a(t) x, where a(t) = 
(1 + z) _1 is the cosmic scale factor, and z the correspond- 
ing red-shift (ao = 1). Thus all geometrical aspects of this 
universal line element are determined up to the function 
a(t) and the arbitrary constant fc, which defines the usual 
open, Einstcin-deSitter and closed universes. The function 
a(t) can be found for each energy content by solving the 
corresponding equations of motion, eg the gravitational 
field equations. 

In this section we consider Einstein's equations: 



Rftv + Ag^ = -8irG (t^ u - - g^T 



(3) 



where T = g^ v T^ v is the trace of the energy-momentum 
tensor; we have included a cosmological constant term to 
keep the equations general at this stage. For an ideal fluid, 
we have: 



pg^v + (p + p) u^u„ 



(4) 



We can now use the field equations and the above 
energy-momentum to find the scale factor a(t) in the met- 
ric: 
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In the fluid approximation, deviations from the mean 
background p are characterized by fluctuations in the den- 
sity and velocity fields. The continuity equation for a non- 
relativistic fluid is (Peebles 1993): 

^- 5(x, t) + V • {[1 + 5(x, t)] v(x, t)} = (7) 
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where 6(x,t) = p{x,r)/~p — 1 is the local density contrast, 
v(x,t) the peculiar velocity (see Eq.[fLl| below), and r the 
conformal time defined by 



dt did 

dr = — — 4=> — = — 

a[t) dt a dr 



The continuity equation [q] can also be written 



0, 



V • v 



(8) 



(9) 



In order to find an equation of motion for the density 
contrast alone we shall resort to the Raychaudhuri equa- 
tion (see eg Wald 1984) 



^ + i 6 2 = -o-y^ + + fl^ u»u v (10) 

where = V^u* 1 , cr^ is the shear tensor, Wy the vorticity 
tensor, and the Ricci tensor, and s the proper time 
parameter; u M is the fluid's 4- velocity, vP = 1, and 



d(i) a; + u(x, i) 



(11) 



It is important to stress that Raychaudhuri's equation, 
Eq. Jl0[ , is purely geometric: it describes the evolution 
in proper time of the dilatation coefficient O of a bundle 
of nearby geodesies. There is no physics in this equation 
until a relationship between i? M „ and the matter contents 
of the universe is specified by means of a set of field equa- 
tions. This makes it very useful for our purposes in this 
paper, as we shall later make reference to a different set of 
field equations. 

If Einstein's field equations, Eq. |3| and [Q, are assumed 
then it is readily verified that 



u ^u y = -AnGp 1 
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2.2. Shear free and matter domination 



In a matter dominated regime (p = 0), p ~ a~ 3 . Equa- 
tion H for the Hubble rate H , can be rewritten using the 
notation: fijw = 8ttGpq/(3Hq), which is the ratio of the 
current matter density to the critical density, Vt^ = k/ Hq 
gives the global curvature, and Ha = A/(3Hq) where A is 
the cosmological constant, so that f2jw +fijt + J^a = 1 ; 



d0 



-6» 2 = -4irGa 2 p6 



(16) 



where 7i(r) = rffln a) /dr. We can now eliminate 6* be- 
tween eqs. || and Jl(| to find the following second order 
differential equation for the density contrast: 



dr 2 
4 



H{r)f T -\n\T)n M {r)5 



(i + sy 



+ -'H 2 {t)Vl m (t)8 2 (17) 



where we have shifted to the rhs all non-linear terms, and 
used the notation 
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Qm + a fife + a 3 £Ia 



(18) 



Equation p7[ reproduces the equation of the spherical 
collapse modeT(SC). In other words, t/ie >9C approximation 
is the exact dynamics when shear is neglected (see Fosalba 
& Gaztanaga 1998a). As one would expect, this yields 
a local evolution, in the sense that the evolved field at a 
point is just given by a local (non-linear) transformation 
of the initial field at the same point, with independence 
of the surroundings. This SC solution yields the exact 
perturbation theory predictions for the cumulants at tree- 
level (leading order with Gaussian initial conditions) and 
it also is an excellent approximation for next to leading 
orders, see below. As mentioned in the introduction, one 
can also use the SC model to predict the value of the crit- 
ical linear overdensity, S c , that will collapse into virialized 
halos. 

2.3. Linear growth 

We next do a perturbative expansion for S. The first 
contribution is the linear theory solution. For this, equa- 
tion [0 clearly simplifies to 



=£ + *w£-§*V>n-M 







(19) 



where Si stands for the "linear" solution. Because the co- 
efficients of the above equation are time dependent only, 
the spatial and temporal part factorise: 



h 2 {z) = h 2 q [ n M (i + z) 3 + n k {i + z) 2 + n A ] (13) 

We can now replace equation into equation |Io| . 
In a matter-dominated regime, and for a shear free, non- 
rotating cosmic fluid we obtain: 



dQ 1 



e 2 



-AnGp + A 



dt ' 3 

On making use of equation |ll|] we can split as 

3a 6 
a a 



e = v M u" 



(14) 



(15) 



so that, taking into consideration the field equations for 
the expansion factor a(t) (Eqs. (|) and pj), equation JL4| 
can be recast in the form 



5 1 (x,t)=5 (x)D(t) 



(20) 



where D is usually called the linear growth factor. Thus 
initial fluctuations, no matter of what size, are amplified 
by the same factor, and the statistical properties of the ini- 
tial field are just linearly scaled. For example, the TV-point 
correlation functions are: 



£/v(n, -.,rN,t) = D N £ N (ri, ..,r Nl 0) 



(21) 



To find the solution to equation fll9| ] it is expedient to 
change the time variable to r\ = ln(a)7so that 



1 d 



dr/ H(t) dr H dt 



(22) 



We then have 
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Fig. 1. — The non-linear density contrast, 8, as a function of the linear one <5; in the spherical collapse. The continuous line shows the 
GR prediction (w = oo, e = 0, 7 = 0), the short-dashed lines correspond to the BD model with u> = 10 and ui = 1 (from left to right). The 
long-dashed line shows the case with a non-standard Hubble rate H 2 = a~ 3 ( 1 + e ) for e = 0.5. 



d 2 D 
d 2 r\ 

where we can write 
H 

BP ~ ' 




n M (r))D = 



Qm + 2/3 e" Qk 
Q M + e" flk + e 3r > fl A 



n 



M 



(23) 



(24) 



(25) 



where SIm, Qk and Qa are just constants (the current value 
at a = 1). 

In the Einstein-deSitter universe (Qk — = 0) we 
have that Qm(v) = 1 an d H / H 2 = —3/2, so the differen- 
tial equation becomes 



d 2 D 1 dD_ 3 D _ 
d^" + 2 _ 2 



(26) 



whose solutions 



D = Ci e" + Cae" 3 / 2 " = d a + C 2 a~ 3 ^ (27) 

reproduce the usual linear growth D ~ a and the decaying 
solutions £) ~ a~ 3 / 2 . 

2.4. Non-linear growth 

The exact (non-perturbative) solution for the SC Eq. p7[ 
for the density contrast in an Einstein-deSitter universe 
admits a well known parametric representation: 

9 Qg- siny) 2 



2(l-cos<p) E 



-(</? — sin tp) 



-1 2/3 



(28) 



for 5; > 0, linear overdensity, and 



9(sinh<^-<^) 2 



2 (coshtp- l) 3 



-(sinh<^ — tp) 



2/3 



(29) 



for 5; < 0, linear under-density (see Peebles 1993), where 
the parameter tp is just a parametrisation of the time co- 
ordinate. There is also a solution for the Qm 7^ 1 case 
(see Bernardeau 1992, Fosalba & Gaztahaga 1998b). The 
continuous line in Figure [I] illustrates the solution to the 
above equation (the other lines will be explained later). 
Note the singularity at Si ~ 1. 686, w hich corresponds to 
the gravitational collapse (see § |3.1.3 below). 

If we are only interested in the perturbative regime 
(Si — ► 0), which is the relevant one for the description 
of structure formation on large scales, the above solution 
can be expressed directly in terms of the linear density 
contrast, Si, which plays the role of the initial size of the 
spherical fluctuation in Eq. p0| . This way, the evolved den- 
sity contrast in the perturbative regime is given by a local- 
density transformation of the linear density fluctuation, 



5 = m) = E S [SlY 

n—l 



(30) 



Notice that all the non-linear dynamical information in 
the SC model is encoded in the v n coefficients. We can 
now introduce the above power series expansion in Eq. jL7| 
and determine the v n coefficients one by one. Before we 
do this, it is convenient to change again the time variable 
to rj = In (a) as we did in the linear case, Eqf23]: 
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Fig. 2. — The linear growth index a± (continuous line) and non-linear coefficients U2 (short-dashed) and V3 (long-dashed), as a function of 
7 = p/p. Vertical dotted lines correspond to the vacuum, matter and radiation dominated cases 7 = —1,0, 1/3. The horizontal dotted lines 
bracketted the 1/2 and 1/3 regions within 10% error of the matter dominated (7 = 0) case. 




(31) 



We can now use the expansion in Eq. |3Cj] with 5i given 
by the linear growth lactor D = a = e n and compare order 
by order. For the Einstein-deSitter universe they turn out 
to be: 
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34 
21 



^3 



682 
189 



(32) 



and so on (see eg Folsalba & Gaztahaga 1998b for other 
cases) . Once we have these coefficients we can get the evo- 
lution of the non-linear variance and hi gher order moments 
in terms of the initial conditions (see §4.1 below). 



2.5. Equation of state p = jp 

We will now consider a perfect fluid with equation of 
sate p = 7/9. Not all values of 7 make physical sense. Here, 
in the spirit of going beyond the standard paradigm, we 
will ignore these restrictions and assume that 7 can take 
any real constant value, irrespective of other cosmological 
parameters. 

The time-component of the energy conservation equa- 
tions — gives us (for p = jp) both the back- 
ground density behavior 

pa 3(1+7) = const (33) 
and the continuity equation for the density contrast 

— + (l + 7 )(l + £)0 = -7}5(t>.W) (34) 
ar 



where, like before, r is the conformal time, and 9 = V • v. 
This is the generalization of equation |^] for a relativistic 
fluid. Note that an additional (quadratic) term now ap- 
pears in the rhs of . The magnitude of this term is as- 
sessed by resorting to the space-components of the energy 
conservation equations V v T^ v = 0: these are identically 
satisfied when 7 = 0, and they show that v-V<5 oc |v| 2 /c 2 , 
plus higher order contributions. These can be safely ne- 
glected since peculiar velocities are always very small com- 
pared to the speed of light; in fact the approximation 



2 /c 2 



is always made, even in the more standard 



case when 7 = 0. We shall therefore consistently adopt 
the following equation for the density contrast: 



^- + (1+ 7 )(1+S)0 = O 
ar 

Also, Hubble's equation, Eq. 113], now becomes 



H 2 = Hi 



-3(1+7) 



(35) 



(36) 



We can combine equation |35| with the Raychaudhuri 
equation for this case — cf Eqs. ]|10| and pa 



■ e 2 



do 1 

~dt + 3 
to obtain, after some algebra 



4wGp (1 + 3 7 ) + A 



(37) 
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2 + TP g-f(l+7)(l+37)fW 



where we have expediently redefined n(rj) in Eq. to 



A I 



e n(i+37) fl k + e 3r i(i+~<) fl A 



(39) 
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Fig. 3. — The critical value of the linear density contrast <5 C where <5 = oo as a function of 7 = p/p. 



and we can write 



H 
IP 



(40) 



In an Einstein-deSitter universe (fife = Q\ = 0), 0(77) = 1, 
and the linear regime is governed by 

d 2 D I-37 dD 3, 
which has the usual solutions of the form D — a Q , with 



ax = l + 3 7 , a 2 = -3(l+7)/2 



(42) 



Figure g shows these perturbative solutions. The shaded 
region corresponds to the case where linear evolution is 
suppressed, eg a < 0. In this case, as can be seen from 
Eq.|43[|-|^9t, v 2 and ^3 have a very rapid variation. The 
growing mode for 7 > —1/3 is: 



«i = 1 + 3 7 

2 (17 + 48 7 + 277 2 ) 



V-2 



V3 



3 (1+7) (7+15 7 ) 
72 + 540 7 + 324 7 2 
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(6 + I87) (17 + 487 + 277 2 ) 
(1+7) (7+15 7 ) 

x (27+ 144 7 + 189 7 2 )" 1 



24 
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(43) 
(44) 

(45) 
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For 7 < — 1 the dominant linear growth is a 2 and the 
values of v 2 and v% are constant: 
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'3 (I + 7) 



(47) 



l>2 = — 

2 

1/3 = 3 



(48) 
(49) 



For radiation (7 = 1/3) we have that ct\ — 2 which 
reproduces the well known results (see Peebles 1993) and 
v 2 — § and 1/3=3, which are new results as far as we 
know. Note that these values are identical to the case of 
negative pressure, 7 < — 1, the only difference being in the 
linear growth, but for 7 = —7/3 all a, vi and ^3 are iden- 
tical to the radiation case. In the limit of strong pressure 
7 — > 00 we find: v 2 = 6/5 and z/3 = 12/7 . As can be 
seen in Figure and also in the equations above, there 
are poles for v 2 at 7 = —1 and 7 = —7/15. 

Figure ^ shows the corresponding variation in 6 C , de- 
fined as the value of the linear overdensity where t he cq r- 
responding non- linear value becomes infinity (see § |3.1.3 ). 



3. GRAVITATIONAL GROWTH OUTSIDE GR 
3.1. Scalar-Tensor Theories 

Here we investigate how a varying G could change the 
above results. We parameterize the variation of G us- 
ing scalar-tensor theories (STT) of gravity such as Brans- 
Dicke (BD) theory or its extensions. 

To make quantitative predictions we will consider cos- 
mic evolution in STTs, where G is derived from a scalar 
field (f> which is characterized by a function to = oj{4>) de- 
termining the strength of the coupling between the scalar 
field and gravity. In the simplest BD models, oj is just a 
constant and G ~ <f)~ l — see below. However if u> varies 
then it can change with cosmic time, so that oj — lo(z). 
The structure of the solutions to BD equations is quite 
rich and depends crucially on the coupling function ui{4>) 
(see Barrow & Parsons 1996). 

Here we shall be considering the standard BD model 
with constant uj\ the field equations are (see eg Weinberg 
1972): 
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n _ (rp 1 + ^ -r 



3 + 2w 
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8tt 



3 + 2w 



T , (T = fl ""T M1/ ) 



(50) 
(51) 



The Hubble rate H for a homogeneous and isotropic 
background universe can be easily obtained from the 
above; 



8np k A u)(f> 2 6 
-^ + ^ + 3 + 6^- H ^ (52) 



These equations must be complemented with the equa- 
tion of state for the cosmic fluid. In a flat, matter domi- 
nated universe (p = 0), an exact solution to the problem 
can be found: 



G=p-^4>- 1 = G,(l + zfl^ (53) 



2uj 



and 



a(t) = ( t /f )(2"+2)/(3"+4) 



(54) 

This solution for the flat universe is recovered in a gen- 
eral case in the limit t — ► oo, and also arises as an exact 
solution of Newtonian gravity with a power law G oc t n 
(Barrow 1996). For non-flat models, a(t) is not a simple 
power law and the solutions get far more complicated. To 
illustrate the effects of a non-flat cosmology we will con- 
sider general solutions that can be parametrized as Eq. 

but which are not simple power-laws in a(t). In this 
case, it is easy to check that the new Hubble law given by 
Eq. [R2J becomes 



H 2 



H 2 



n M (i + z ) 3+1 ^ 1+ ^ + h k (i + z f + Ci A 



where 0^,0^ and £l A follow the usual relation Qm + &k + 
fl A = 1, and are related to the familiar local ratios (z — > 0: 
U M = 8^G po/(3i? 2 ), n k = k/H§ and fi A = A/(3ff 2 )) 
by 
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(56) 



Thus the GR limit is recovered as to — * oo. 

We now investigate the density fluctuations in the above 
theory. Like in section II, we shall make use of the conti- 
nuity equation |^] in combination with the Raychaudhuri 
equation fl(i|| . As mentioned above, cf section 2.1, both 
of these are still valid within the context of BD theory: 



it is only needed to replace the Ricci tensor in the rhs of 
Eq. 1 3-0 1 according to BD's field equations, Eq. |5C|] . Con- 
sidering again a non-rotating, shear-free cosmic fluid, we 
find: 



d& 
~dt 



e 2 



4 + 2uj Airp f I 



u> 3p 



2 + lo p 



(57) 



We shall still make use of a g ravitational "constant" 
parametrized as in equation [p3|| above; this is justified 
insofar as the characteristic length for the variation of <fi 
is typically much greater than that of the density fluctu- 
ations in a matter dominated universe — see eg (Nariai 
1969). In this approximation, the above equation gives 



d0 „„ , 



4 + 2uj 4:Tra 2 p5 
3 + 2uj 6 



(58) 



where r is again the conformal time parameter, dr 
a dt, and 9 is defined in equations |)| and Jig ]. Remark- 
ably, this equation is very similar to the GR equation 
we only need to replace in it the gravitational constant G 
by its expression as a multiple of the varying scalar field <fi 
given in equation ]55( ] . Combining |58j with the continuity 
equation B we immediately find 



d 2 5 



dS 



dr 2 dr 3 (1 + 8) \ <h I ■', + ±j 



dS 



A + 2u 4ira 2 p6 



(59) 



Like in section II, we change the independent variable 
in p9| to rj = In a, whereby we obtain 




(1 



sy 1 ( 



MY 
d^J 



2uj H 2 (j) 



(60) 



Using equation |5^] to calculate H, and assuming fur- 
ther that Ct k = — 0, we finally get 



d 2 S 1 



dS 1 (2 + w)(4 + 3w) 



d 2 r) 2 1 + u di) 2 (1 + lo) 2 

2 



4 1 



3 1 



dS 
6 V drj 



1 (2 + ^)(4 + 3u;) 2 

2 (l + to) 2 



(61) 



We next examine the solutions to this equation. 

3.1.1. Linear growth 
Let us call D(j]) the solution to the linearized version of 



equation |61|, i.e. 



d 2 D 1 



dD 1 (2 + w)(4 + 3w) 



d 2 rj 2 1+uj drj 2 {l + tof 



D = (62) 



Again the solutions are given by the roots ct\ and a\ of 
the corresponding characteristic functions: 
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Fig. 4. — The linear growth indices «i (continuous line) and &2 (dashed line), defined by the solution D ■■ 
the BD parameter ui for a time varying gravitational constant G = Go a~ . 



Ci a" + C2 a" as a function of 



D = d a ai + C 2 a° 



with 



■ uO 



Oil 



a 2 



1+LO 

-4 -3a; 
2 + 2u 



l - + o(±- 2 

L0 Z 



uO 



2 lo 



LO' 



(63) 

(64) 
(65) 



which reproduces the usual linear growth D ~ a and 
D ~ a -3 / 2 in the limit lo — > 00. Note that ai corre- 
sponds to the growing mode only for large values of 
but the situation is more complicated when lo is not large. 

Figure ^ shows the values of ai and «2 as functions of 
lo. The effective G in BD decreases as the Universe ex- 
pands if — 1 < lo < 00, and the expansion factor a(t) stops 
for lo = — 1; the growing mode in this regime is controlled 
by ai, since this is the positive root. The growing mode 
for —4/3 < lo < —1 is 012, but the universe shrinks to an 
eventual collapse in this regime (see Eq. J54|). Between 
—2 < lo < —4/3 the Universe expands again, but there 
are no growing modes, as can be seen in Figure |] (both ct\ 
and a2 are negative). For lo < — 2 the expansion factor 
grows with time and a\ becomes the growing mode again. 
Notice that in this regime of lo < —2, a% < 1, so that it 
is slower than for lo > 0. As we will show below this is 
compensated in part by a stronger non-linear growth. 

3.1.2. Non-linear growth 

In the non-linear case we consider the full version of 
equation We can now proceed as before, using the 

expansion in Eq. ]30| with Si given by the linear growth 
factor D = a ai — e ain , and compare order by order. We 
find 



V-2 



Mlo + 56 
21cj + 36 



34 
21 



8 1 
119 u> 



O 



(66) 



2(944+ 1136 lo + 341 lo 2 ) 
3(12 + 7w)(16 + 9w) 



682 
189 



3452 



21483 lo 



(67) 



Note how for positive lo non-linear effects tend to com- 
pensate the increase in linear effects, cf Figure |], whereas 
for lo < —4/3, the linear effects are reduced (a < 1) while 
non-linearities get larger. 

Figure H shows the variation in 1/2 as a function of lo 
using Eq.[|66[. Negative values of lo produce almost sym- 
metrical variations in the opposite direction when |u;| is 
large. For small lo there is a pole at lo = —12/7 where V2 
diverges. But note that there is no growing linear mode in 
this case, which means that fluctuations are rapidly sup- 
pressed. 



3.1.3. Strongly non-linear regime 

Figure |l| shows the fully non-linear solution for the over- 
density S as a function of the linear one Si. The contin- 
uous line shows the standard solution to Eq.jl7| as given 
in Eq.[p8[-[p9[. As can be seen in the Figure, there is a 
critical value of 61 = 3/2(3?r/2) 3 / 2 ~ 1.6865 where the 
non-linear fluctuations become infinite. This corresponds 
to the point where the spherical collapse occurs (see Pee- 
bles 1993). Thus an initial fluctuations So will collapse 
after evolving a time t, such that the growth factor is 
D{t) = S c /5o- For the standard GR, flat and matter dom- 
inated case, this time would correspond to a formation 
red-shift: Zf = 6q/S c — 1 (if we use a — 1 today). For the 
BD case both S c and D(t) are different, so that formation 
times Zf will be correspondingly different (see Eq.flsll). 
The short-dashed lines in Figure [P correspond to the same 
exact solution in the BD model with lo = 10 and lo = 1. 
Right panel in Figure || illustrates how 5 C changes in the 
BD model as a function of lo. 
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Fig. 5. — Left panel: dashed lines show U2 (top) and 1/3 (bottom) as a function of u> for a time varying gravitational constant 
G = Go a _1 /' 1+ ^\ The GR results, G = Go, (continuous horizontal lines) are bracket by 10% errors (dotted lines). Right panel: the 
critical value of the linear density contrast <5 C where 8 = 00 as a function of u). 



3.2. Gravitational Growth with H 2 



-3(l + £) 



Consider now the flat case with Q,k = £1a = 0. To 
account for a simple variation on the standard Einstein's 
field equations we will consider the case where fluctuations 
grow according to the matter dominated case (ie 7 = 0) 
but the background evolves in a different way. We will as- 
sume that the Hubble rate goes like H 2 ~ a- 3 ( 1+e ) rather 
than H 2 ~ a~ 3 . It might be possible to find some moti- 
vation for this model, but this is beyond the scope of this 
work. Here we just want to introduce some parametric 
variations around the standard field equations to see how 
things might change. In this case we have: 



d 2 5 l-3e dS 3 4 
d 2 7] + 2 dq ~ 2 ~ 3T 



(68) 



The solutions for the linear growth factor index and the 
non-linear coefficient V2 are 



«i = 



-1 + 3e + V25 - 6e + 9e 2 



131 - 30e 



4 

45e 2 



(1 - 3e) V25 - 6e + 9e 2 



84 - 18e + 27e 2 



(69) 
(70) 



These solutions as a function of e are illustrated in Fig- 
ure ||, which also shows v^. As can be seen in the Figure, 
the higher the linear growth index ol\ the lower the non- 
linear coefficients. 

Right panel in Figure ^ shows the corresponding varia- 
tion in S c . 

4. OBSERVATIONAL CONSEQUENCES 

We will focus here on Gaussian initial conditions. That 
is, our initial field for structure formation is a spatial real- 
ization of a (three-dimensional) Gaussian distribution with 



a given power spectrum shape, and a very small initial am- 
plitude. As we are interested in the gravitational regime 
alone, this field will be smoothed over a large enough 
scale, corresponding to the distance beyond which non- 
gravitational forces (eg hydrodynamics) can be neglected. 
Thus at each point the overdensity 6(x) grows according 
to gravity, which in the shear free approximation is just a 
local dynamics: the spherical collapse (eg Eq.|l7||). 

4.1. Cumulants 
Consider the J — order moments of the fluctuating field: 



TO J 



(71) 



Here the expectation values ( ... ) correspond to an av- 
erage over realizations of the initial field. On comparing 
with observations we assume the fair sample hypothesis 
(§30 Peebles 1980), by which we can commute spatial in- 
tegrals with expectation values. Thus, in practice (■••) 
is the average over positions in the survey area. In this 
notation the variance is defined as: 



Var(5) — a 2 = n%2 — m\ 



(72) 



More generally, we introduce the connected moments 
£j, which carry statistical information independent of 
the lower order moments, and are formally denoted by 
a bracket with subscript c: 



(73) 



The connected moments are also called cumulants, re- 
duced moments or irreducible moments. They are defined 
by just subtracting the lower order contributions: 
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Fig. 6. — Left panel: The linear growth index a\ (continuous line) and non-linear coefficients V2 (short-dashed) and 1/3 (long-dashed), as a 
function of e, which parameterizes a non-standard Hubble rate H 2 ~ a -3 ( 1+e ). Vertical dotted line corresponds to the standard Hubble law 
(e = 0). The horizontal dotted lines bracket the V2 and ^3 regions within 10% error of the standard (e = 0) case. Right panel: the critical 
value of the linear density contrast (5 C where 8 = 00 as a function of e. 



£j = mi = 



-2 



£ 2 = er = m 2 - £1 = m 2 
£3 = m 3 - 3£ 2 £i - ?! = m 3 



(74) 



£4 = m 4 - 4f 3 f 1 - 3£ 2 - 6£ 2 £i - £ x = m 4 - 3m 2 
and so on. It is useful to introduce the hierarchical ratios: 

h 



Sj 



(75) 



which are also called normalized one-point cumulants or 

reduced cumulants. We shall use the term skewness, for 
2 3 

S3 = £3/^2 an d kurtosis, for S4 = £4/^2- 

4.1.1. Linear Theory 



As mentioned in §2.3, initial fluctuations, 5q, no matter 
of what amplitude, grow all by the same factor, D; thus 
the statistical properties of the initial field are just linearly 
scaled in the final (linear) field, Si: 

(6f) e =D J {5{) c (76) 

Consider for example the linear rms fluctuations 07 or 
its variance of. In the linear regime we have: 



of = (S 2 (t)) = (D(t-t ) 2 S 2 )=D(t-t ) 2 o 2 (77) 

where errj refers to some initial reference time to. To give 
an idea of this effect, consider the growth of fluctuations 
since matter domination, when the universe was about 
1100 times smaller. In General Relativity (GR) in the mat- 
ter dominated Einstein-deSitter universe, a would grow by 
a factor D ~ 1100. While, if we take w ~ 10 in the DB 



theory, eg Eq.[g4j, we have that fluctuations increase in- 
stead by a factor D ~ 2079, which is about 1.9 times larger 
in a, so the variance nowadays would be about 3.6 times 
larger if we fixed it around the COBE variance of Cosmic 
Microwave Background (CBM) temperature fluctuations. 
For us ~ 100, the variance would only be 14% larger than 
in GR. This latter result is small, but it could be relevant 
for future precision measurements (eg MAP or PLANCK 
satellites to map CMB and 2DF or SLOAN DIGITAL SKY 
galaxy surveys). Similar considerations can be made for 
the values of a with a different cosmic equation of state, 
eg Eq.[[l3| or a different Hubble law, Eq.|6£|. In general 
we can write that a small change in a would produce a 
relative change in the linear rms of 



Ac- 
er 



ln(l + z)Aa 



(78) 



Thus, a change of only 1% in the absolute value of the 
equation of state 7, would produce a relative change of 
20% in a between recombination z ~ 1100 and now, cf 
Eq.@. 

The hierarchical ratios (see Eq.(75[) will scale as, Sj = 
Sj(0)/D J ~ 2 , where Sj(0) are the initial ratios. This im- 
plies that the linear growth erases the initial skewness and 
kurtosis, so that Sj — * 0, as time evolves (and D — * 00). 
Note that if we want to do a meaningful calculation of 
these ratios or the cumulants, in general we might need 
to consider more terms in the perturbative series, Eq. |p(|. 
For Gaussian initial conditions Sj(O) = 0, and we need to 
consider higher order terms in the perturbation series to 
find the leading order prediction. 

4.1.2. Weakly non-linear 

The next to leading order solutions for the cumulants of 
the evolved field given the expansion Eq. fl3Cl] , can be easily 
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found by just taking expectation values of different pow- 
ers of S (see eg Fosalba & Gaztafiaga 1998a). For leading 
order Gaussian initial conditions we have 



S 3 = 3v 2 + 0(af) 

Si = \v 3 + 12t| + 0(af) 



(79) 



For non-Gaussian initial conditions see Fry & Scherrer 
(1994) Chodorowski & Bouchet (1996), Gaztafiaga & Ma- 
honen (1996), Gaztafiaga & Fosalba (1998). 

If we use for v 2 the solution in Eq.^], eg v 2 — 34/21, 
the skewness yields S3 = 3v 2 = 34/7, which repro- 
duces the exact perturbation theory (PT) result by Pee- 
bles (1980) in the matter dominated Einstein-deSitter uni- 
verse. Thus the shear-free or SC model gives the exact 
leading order result for the skewness. This is also true 
for higher orders (see Bernardeau 1992 and Fosalba & 
Gaztafiaga 1998a) and for other cosmologies (eg Bouchet 
et al. 1992, Bernardeau 1994a, Fosalba & Gaztafiaga 
1998b, Kamionkowski & Buchalter 1999). For smoothed 
fields, the exact leading order results are slightly different: 



S3 
S 4 



34 

y 

60712 62 



•7i 



1323 



■7i 



■7i 



(80) 
(81) 



where 71 is the logarithmic slope of the smoothed variance 
(see Juszkiewicz et al. 1993, Bernardeau 1994a, 1994b). 
These can also be reproduced in the shear-free approx- 
imation as shown by Gaztafiaga & Fosalba (1998); this 
results in a smoothing correction: 



V2 



V2 



21 

3 



v 3 = -(-271 



+ 71 + 671 v 2 + 4^) 



(82) 



and replacing v 2 and 1/3 by V2 and P3 in Eq. |^] ( see Fos- 
alba & Gaztafiaga 1998a for more details). There are also 
corrections to the above expressions when measurements 
are taken in red-shift space (eg Hivon et al 1995, Scoc- 
cimarro, Couchman and Frieman 1999). Next to leading 
order terms have been estimated by Scoccimarro & Frie- 
man (1996) (see also Fosalba & Gaztafiaga 1998a, b). 

The smoothed values of 63 and 64 can be measured as 
traced by the large scale galaxy distribution (eg Bouchet et 
al. 1993, Gaztafiaga 1992, 1994, Szapudi el at 1995, Hui 
& Gaztafiaga 1999 and references therein), weak-lensing 
(Bernardeau, Van Waerbeke & Mellier 1997, Gaztafiaga 
& Bernardeau 1998, Hui 1999) or the Ly-alpha QSO ab- 
sorptions (Gaztafiaga & Croft 1999). These measurements 
of the skewness and kurtosis can be translated into estima- 
tions of v 2 and V3 which can be used to place constraints 
on 7, ui or e using Eq.Q], Q and For small values 
of these parameters the relationship is linear, so the un- 
certainties in S3 and S4 would directly translate into the 
corresponding uncertainties in 7, lu or e. 

The expressions above apply to unbiased tracers of the 
density field; since galaxies of different morphologies are 
known to have different clustering properties, at least some 



galaxy species must be biased tracers of the mass. As an 
example, suppose the probability of forming a luminous 
galaxy depends only on the underlying mean density field 
in its immediate vicinity. Under this simplifying assump- 
tion, the relation between the galaxy density field S ga i(x) 
and the mass density field 6(x) can be written as 



(83) 



where b n are the bias parameters. Thus, biasing and grav- 
ity could produce comparable non-linear effects. To lead- 
ing order in £ 2 , this local bias scheme implies ^ 9 2 = b\£ 2 , 
and (see Fry & Gaztafiaga 1993) 



gal 



S 



ggal 



^ + 3 ^ 



1 

b 2 S 3 



b\ ^ b\ 



63 bl 
44 + 12^f 



(84) 



Gaztafiaga & Frieman (1994) have used the comparison 
of S3 and S4 in PT with the corresponding values mea- 
sured APM Galaxy Survey (Maddox et al. 1990), to infer 
that 61 ~ 1, i>2 — and 63 ~ 0, but the results are de- 
generate due to the relative scale-independence of Sjy and 
the increasing number of biasing parameters. One could 
break this degeneracy by using the configuration depen- 
dence of the projected 3-point function, 53(a), as proposed 
by Frieman h Gaztafiaga (1994), Fry (1994), Matarrese, 
Verde & Heavens (1997), Scoccimarro et al. (1998). As 
shown in Frieman & Gaztafiaga (1999), the configuration 
dependence of q 3 (a) on large scales in the APM catalog is 
quite close to that expected in perturbation theory (see 
Fry 1984, Scoccimarro et al. 1998, Buchalter, Jaffe & 
Kamionkowski 2000), suggesting again that 61 is of order 
unity (and b 2 — 0) for these galaxies. These agreement 
indicates that large-scale structure is driven by non-linear 
gravitational instability and that APM galaxies are rela- 
tively unbiased tracers of the mass on these large scales. 

The values of S3 and S4 in the APM are measured 
to agree with the standard matter dominated Einstein- 
deSitter universe within about 10% — 20% (see Gaztafiaga 
1994; Gaztafiaga & Frieman 1994; Baugh, Gaztafiaga & 
Efstathiou 1995; Gaztafiaga 1995, Hui & Gaztafiaga 1999), 
also in agreement with the shape information in the 3- 
point function (see Frieman & Gaztafiaga 1999). For 
example, using the projected APM catalogue Gaztafiaga 
1994 (Table 3) finds an average of S 3 = 3.2 ± 0.2 and 
S 4 ± 20.6 ± 2.6 scales between 7 and 30 h~ l Mpc. For an 
average APM slope of 71 ~ 1.7, these values are in agree- 
ment with the PT predictions in Eq.|yJ yield: S3 ~ 3.1 
and S 4 ~18. 

The 1-sigma error bar of ~ 10% on large scales quoted 
by Gaztafiaga 1994 is mostly statistical (sampling error). 
Other systematics effects due to biasing, projection, or 
large scale errors in the building of the APM catalogue 
could be of the same order (see Frieman & Gaztafiaga 
1999 and Hui & Gaztafiaga 1999). Thus given the cur- 
rent uncertainties it would be conservative to take a 20% 
error bar. Unfortunately, with such large error bars we 
can not constraint much the values of 7 uj or e. Stronger 
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constraints can be found if we take the more optimistic 1 
sigma 10% error bars in the measurements of S3 and S. 
This case is shown as horizontal dotted lines in Figures 
and |(]. From v 2 the 10% uncertainty translates into 



I 



0.2 < 7 < 0.4 
-2.4 > lj > -1.0 
-0.9 < e < 0.9 



(85) 



Note that this is still of marginal interest. For example, 
the constraints on 7 include the possibility of a radiation 
(7 = 1/3), matter (7 = 0) or negative pressure 7 < 0. 
From 1/3 we can obtain stronger constraints from a 10% 
error (but obviously systematic effects could be larger for 
higher order cumulants): 



-0.1 < 7 < 0.15 
-3.4 > cu > -0.2 
-0.35 < e < 0.35 



(86) 



These bounds are more interesting. It is clear that 
forthcoming surveys (such as the SLOAN Digital Sky Sur- 
vey) will dramatically improve this situation (for errors on 
statistics see Szapudi, Colombi and Bernardeau 1999, and 
references therein). 

Note that the above results are independent of the nor- 
malization of fluctuations. 

4.2. Collapsed objects 

Press & Schechter (1974) formalism and its extensions 
(eg Bond et al. 1991; Lacey & Cole 1993) predict the 
evolution of the mass function of halos and also their clus- 
tering properties. Comparison with N-body simulations 
show a very good agreement of these prescriptions for a 
wide range of statistical properties (eg see Lacey & Cole 
1994 and references therein). For example, the comoving 
number density of collapsed objects (halos or clusters) of 
mass M is 



n(M)dM = -^ 



din a 
din A/ 



exp 



JL 

2a 2 



pdM 
M 2 



(87) 

where a = o~{R) is the current linear rms fluctuation at 
the scale R corresponding to the mass M — 4/3nR 3 p, and 
p is the mean background. The value of S c corresponds to 
the value of the linear overdensity at the time of collapse. 
The collapsing structure virializes when the (non-linear) 
overdensity becomes very large (S > 100). The actual 
definition is not very important, as once d > 100, the non- 
linear collapse is quite rapid, as can be seen in the plots 
of Figure [l| and the corresponding value of Si does not 
change much. Here we will take S c to be the critical value 
where S — > 00; other prescriptions (eg the value of Si corre- 
sponding S ~ 178) yield similar results. For the standard 
Einstein-de Sitter case we have S c ~ 1.686. Note that the 
above abundance depends on the ratio 



V — 





The time of collapse or formation is just given by the 
ratio of S c to the linear overdensity Si today 



Si 



l/c 



1 



(89) 



so that an object which has Si = S c now, has a forma- 
tion red-shift Zf = 0, while a fluctuation 4 times larger 
collapses at Zf = 3 if a = 1 or at Zf = 1 if a = 2. 

Non-standard parametrisation of the spherical collapse 
considered in the previous sections can change the above 
formalism in two ways. If we label objects by its initial 
overdensity So then the corresponding Si today is 



Si = S a° 



(90) 



So a different value of a, from the standard GR result 
(Aa = a — acn), as shown in Figures 0, ||, and [| will 
produce a different amp litude of linear fluctuations today. 
Moreover, as shown in § 3.1.3] and Figures [| [|, and |^, the 
solution to the spherical collapse equation produces dif- 
ferent values of S c , and therefore different mass functions 
and formation times. Finally, for a directly measurable 
quantity, such as the surface density of objects, typically 
one needs the volume element, which is also a function of 
the cosmology. 

For example if fluctuations are normalized at a given 
red-shift, z n , then the change in Si today will be 



^ = Aq log(l + z„) 
Si 



(91) 



For recombination, eg COBE normalization, we have 
z n ~ 1100, and 



ASi 
Si 



3 Ac 



(92) 



In the case of the BD theory, we can see in Figure [| 
that for u) > — 1, Aa = a — ctQR > 0, which means that 
ASi > 0. This makes sense as the linear growth is faster 
and, for fixed initial fluctuations, the final linear overden- 
sity will be larger. As shown in right panel of Figure ||, 
S c will also be larger. Thus in this case the effects tend to 
compensate each other. This is true for both the formation 
red-shift 27 or for v in Eq.pq] above. For the formation 
red-shift Zf we have 



A*/ 
l + z f 



1 



Si 



AS C 

~~s7 



(93) 



which is only valid for small changes. In the BD example 
given above with to — 10 (and COBE normalization) we 
have that AS1/S1 ~ 0.9 while AS C /S C — 0.01, so the net ef- 
fect is still quite large. In this case, a formation red-shift of 
Zf — 1 will change to Zf = 1.39. Thus, a positive finite w 
(which corresponds to a larger G at high red-shifts) tends 
to produce larger (earlier) formation red-shifts and higher 
densities (or larger abundances) at a given red-shift, than 
the standard model. This goes in the direction of some 
recent observations (eg see Bahcall & Fan 1998; Robinson, 
Gawiser & Silk 1998, Willick 1999), which seem to need 
larger abundances that expected in some standard cos- 
mologies. This interpretation is degenerate with respect 
to initial conditions and cosmological parameters. 
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for three cosmologies oCDM {M th = 1.9, f2 m = 0.3, Q A = 0, h = 0.65, T = 0.25, erg = 1.0); ACDM (M th = 2.2, O m = 0.3, ft A = 0.7, h = 
0.65, T = 0.25, as = 1.0) and tCDM (M th = 1.3, f2 m = 1.0; Ha = 0, h = 0.5, T = 0.25, erg = 0.56) derived using the Press-Schechter 
prescription. The lower (upper) dashed lines correspond to the a Brans-Dicke Cosmology with lu = 100 {lo = 25) normalized to COBE with 
the tCDM model. . 



Figure [7| illustrates the large differences in the cluster 
counts that can be seen between different cosmological 
models at z > 1 (sec Holder et al. 1999 for details). De- 
viations from General Relativity in the BD models with 
lo = 100 and lo = 25 can be noticed even at low redshift, 
when models are normalized to CMB fluctuations. 

A similar trend is found for the case of Hubble rate 
H 2 = a~ 3 ( 1+c ) parametrisation. A change of |Ae| ~ 0.3 
(allowed by the bounds in Eq.[^6|), when normalized to 
COBE, also produces |A<5//<5;| ~ 0.9 and a smaller effect 
on A5 C /S C . This translates into a similar change (of several 
tens to hundreds of percent) in zj. Earlier (later) forma- 
tion times and larger (smaller) abundances are found for 
e > (e < 0). 

The change in the equation of state p = jp could pro- 
duce comparable effects. The allowed values in Eq.[jS(Sf 
of IA7I ~ 0.1 translate into |A5;/i5/| ~ 0.3, which results 
in similar changes for Zf in either direction, with earlier 
formation for 7 > 0. 

If the normalization is not fixed, ie we do not quite know 
what is the value of the initial fluctuation that gave rise 
to an object we see today (eg a cluster), then all the rela- 
tive change in the formation or abundance comes through 
S c , which tends to produce smaller (later) formation red- 
shifts (S c is larger than the standard GR value) and lower 
densities (or smaller abundances) at a given red-shift. 

5. DISCUSSION AND CONCLUSIONS 

We have reconsidered the problem of non-linear struc- 
ture formation in two different contexts that relate to ob- 
servations: 1-point cumulants of large scale density fluc- 
tuations and the epoch of formation and abundance of 
structures using the Press & Schechter (1974) formalism. 
We have use the the shear-free or spherical collapse (SC) 
model, which is very good approximation for the above ap- 



plications. We have addressed the question of how differ- 
ent are the predictions when using a non-standard theory 
of Gravity, such as BD model, or non-standard cosmologi- 
cal model (eg a different equation of state or Hubble law). 
Note that these are slight variations on the standard theme 
in the sense that they preserved the main ingredients of 
GR, such as the covariance and the geometrical aspects 
of the theory, including the same metric, with only slight 
changes in the field equations. 

We have also presented some preliminary bounds on 7, 
lo and e from observations of the skewness and kurtosis in 
the APM Galaxy Survey, eg Eq.||-@. These bounds 
are optimistic given the current data, but the situation is 
going to change rapidly, and one can hope to find much 
better bounds form upcoming data (such as 2DF or SDSS 
projects). In terms of the equation of state the bounds 
in Eq. |p6| would indicate that our Universe is neither ra- 
diation (7 = 1/3) or vacuum dominated (7 = —1), but 
somewhere in between (eg matter dominated). In terms of 
the Gravitational constant, the bounds on lo from Eq.pof 
would say that G has not changed by more than ~ 5% 
from z ~ 1.15, or by distances of ~ 400/i -1 Mpc. Clus- 
tering at higher red-shift would probe much larger scales 
and times. In terms of e the bounds Eq.[|(|, would say 
that the Hubble law does not differ by more than 7% from 
the standard result (assumed here to be e = 0). We have 
also shown how halo and cluster abundances and forma- 
tion times could change in these non-standard cases. The 
above bounds 0117, w and e from observations of the skew- 
ness and kurtosis in the APM still allow significant changes 
(of several tens to hundreds of percent) on format ion red- 
shifts zj and the corresponding abundances (see §4.2). 

In the context of BD models the limits we find for lo are 
less restrictive than the solar system limits lo > 100. How- 
ever, BD models allow 00 = uo((j)) so that lo can increase 
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with cosmic time, to = w(z), in such a way that it could 
approach the general relativity predictions (ui — > oo) at 
present time and still give significant deviations at earlier 
cosmological times. It is important to recall that our the- 
ory of gravity has only be tested on stellar distances (a.u.) 
while we want to use it on cosmological scales (Mpc). Our 
working example shows, for the first time, how non-linear 
effects are changed in such a model and sets the frame- 
work to study non-linear effects of more complicated (or 
realistic) Scalar- Tensor theories of gravity. 

It is straightforward to combine several of the changes 
proposed here to explore more general situations. One 
could for example parameterize theories in the (7, u>) 
plane, eg different equations of state with different BD 
parameters, or consider the whole (7, u, Om, ^a) space. 
One could also consider a different equation of state for 
the A— component, as in quintessence cosmologies (Cald- 
well, Dave, Steinhardt 1998), such models have already 
been used to predict cluster abundances within the "stan- 
dard" cosmology (see Haiman, Mohr, Holder 2000 and ref- 
erences therein). This would obviously allow for a wider 
set of possible solutions and degeneracies. One should also 
consider other observational consequences of these vari- 
ations, in particular relating to BD theory, such as the 
age of the Universe, the effects on CMB (eg see Chen & 
Kamionkowski 1999), radiation-matter transition (Liddle, 
Mazumdar & Barrow 1998), or the constraints from nu- 
cleosynthesis (Santiago et al. 1997). These considerations 
could rule out some aspects of the proposed variations on 
the standard model, or might require more elaborate so- 
lutions (eg uj = uj(4>) which implies lu — uj(z)). But even 
if this were the case, we still have learn a few new things 
about how structure formation depends on the underly- 
ing theory of Gravity, which is a first step towards further 
analysis of these issues. 

Throughout this paper we have assumed Gaussian ini- 
tial conditions and no biasing. Both biasing (eg Fry & 



Gaztahaga 1993) and non-Gaussianities in the initial con- 
ditions (Gaztahaga & Fosalba 1998) would provide an ad- 
ditional source of degeneracy as they might produce sim- 
ilar effects as the non-standard variations presented here. 
This is the case for example when we have non-zero initial 
skewness or kurtosis, which could produced quite differ- 
ent values of S3 and S4 (eg see Gaztahaga & Mahonen 
1996; Peebles 1999a,b; White 1999; Scoccimarro 2000), 
and therefore to the inferred values of v-^ and V3. Biasing 
can have a very similar effect (eg see Mo, Jing & White 
1997). One would also expect some level of degeneracy 
with biasing and initial conditions for cluster abundances 
or formation times (see Robinson, Gawiser & Silk 1998, 
Willick 1999). 

Rather than proposing an alternative theory of gravity 
or cosmological model, the aim of this paper was to show 
that some small deviations from the current paradigm have 
significant and measurable consequences for non-linear 
structure formation. This could eventually help explaining 
some of the current puzzles confronting the theory, such as 
the need of non-baryonic dark matter. Alternatively, cur- 
rent and upcoming observations of non-linear clustering 
and mass functions can be used to explore our assump- 
tions and place limits on the theory of gravity at large 
(> lh^ 1 Mpc) scales. This provides an interesting test for 
gravity as the driving force for structure formation and 
for our knowledge of the cosmological equation of state. A 
more comprehensive comparison with particular scenarios 
is left for future work. 
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